Key role of hybridization between actinide 5/ and oxygen 2p orbitals 
for electronic structure of actinide dioxides 



m 
o 

(N 
X> 

IX, 



Yu Hasegawa 1 , Takahiro Maehira 2 , and Takashi Hotta 1,3 
1 Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan 
2 Faculty of Science, University of The Ryukyus, Nishihara, Okinawa 903-0213, Japan 
3 Advanced Science Research Center, Japan Atomic Energy Agency, Tokai, Ibaraki 319-1195, Japan 

(Dated: February 13, 2013) 

In order to promote our understanding on electronic structure of actinide dioxides, we construct 
a tight-binding model composed of actinide 5/ and oxygen 2p electrons, which is called f-p model. 
After the diagonalization of the f-p model, we compare the eigenenergies in the first Brillouin zone 
with the results of relativistic band-structure calculations. Here we emphasize a key role of f-p 
hybridization in order to understand the electronic structure of actinide dioxides. In particular, it is 
found that the position of energy levels of I\ and states determined from crystalline electric field 
potentials depends on the f-p hybridization. We clarify the condition on the f-p hybridization to 
explain the electronic structure which is consistent with the local crystalline electric field state. We 
briefly discuss the region of the absolute values of the Slater-Koster integrals for f-p hybridization 
concerning the appearance of octupole ordering in Np02. 

PACS numbers: 71.27.+a, 31.15.aq, 71.15.-m,71.70.Ch 
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I. INTRODUCTION 

Actinide dioxides form a group of important materials 
from technological viewpoints of a nuclear reactor fuel 
and a heterogeneous catalyst. On the other hand, this 
material group has been actively investigated also from a 
viewpoint of basic science because of its high symmetry 
of the fluorite structure of the space group Fm?>m~~— 
In the circumstance of such high symmetry of crystal 
structure, it is possible to observe peculiar ordering of 
multipole higher than dipole, when we change the kind 
of actinide ions. Among several magnetic properties of 
actinide dioxides, a mysterious low-temperature ordered 
phase of NpC>2 has attracted continuous attention in the 
research field of condensed matter physics. 

The phase transition in NpC>2 has been confirmed 
in 1953 from the observation of a peak in the specific 
heat around 25 K 4 Due to the behavior of magnetic 
susceptibility^ it has been first considered that the anti- 
ferromagnetic order occurs. Unfortunately, no dipole mo- 
ments have been observed in the low-temperature phase 
and a puzzling situation has continued. Neutron scat- 
tering experiments have revealed that the ground state 
of the crystalline electric field (CEF) potential is Tg, 7 
which carries multipole moments. Then, from several 
phenomenological works on the ordered phase, a key role 
of octupole degree of freedom has been focused^— In 
fact, the octupole ordering has been strongly suggested 
by 17 0-NMR experiment^ 3 - and also by inelastic neutron 
scattering studyJ^ As for the microscopic origin of such 
higher-order multipole ordering, it has been shown that 
octupole order is stabilized by the orbital-dependent su- 
perexchange interaction, obtained by the second-order 
perturbation of 5/ electron hopping in the Tg degenerate 
Hubbard model on a fee lattice.— — Recently, signifi- 
cant contribution of dotriacontapole moment has been 
also pointed out i 18 i 19 



Since the multipole moments originate from 5/ elec- 
trons, it seems to be natural to consider the Hubbard- 
like model of 5/ electrons. However, from the crystal 
structure of actinide dioxides, it is also important to in- 
clude explicitly 2p electrons, since actinide ion is sur- 
rounded by eight oxygens and the main hopping process 
between nearest neighbor sites should occur from the /- 
p hybridization. In this sense, f-p model is more re- 
alistic Hamiltonian for actinide dioxides. In fact, the 
f-p model for actinide dioxides has been analyzed in 
the fourth-order perturbation theory in terms of f-p 
hybridization^ Then, it has been revealed that octupole 
order actually occurs even when we include oxygen 2p 
electrons. However, there has been a peculiar point that 
the octupole phase appears only for the small absolute 
value of (fpTr)/{fpa), where (fper) and (/fwr) are Slater- 
Koster integrals between / and p orbitals. The reason of 
the sensitivity of the octupole ordered phase concerning 
the f-p hybridization has not been understood yet. 

In order to clarify the role of f-p hybridization for 
the appearance of octupole ordering, Maehira and Hotta 
have performed the band-structure calculations for ac- 
tinide dioxides by a relativistic linear augmented-plane- 
wave method with the exchange-correlation potential in a 
local density approximation^ It has been found that the 
energy bands in the vicinity of the Fermi level are mainly 
due to the hybridization between actinide 5/ and oxygen 
2p electrons. It has been also pointed out that the elec- 
tronic structure at the T point in the first Brillouin zone 
is not consistent with that of the local CEF state. One 
reason for this inconsistency is that the CEF potentials 
are not satisfactorily included in the calculations, but it is 
difficult to control the magnitudes of CEF potential and 
f-p hybridization in the band-structure calculations. It 
is highly requested to reveal the role of f-p hybridization 
for the simultaneous explanation of the octupole ordering 
and the local CEF states. 



In this paper, in order to clarify the roles of hybridiza- 
tion between actinide 5/ and oxygen 2p electrons for the 
electronic structure of actinide dioxides, we analyze the 
tight-binding f-p model in detail. Except for the Slater- 
Koster integrals of (fpir) and (fpcr), we determine the 
parameters in the model from the comparison with ex- 
perimental results and band-structure calculations. In 
order to reproduce the result of the relativistic band- 
structure calculations and obtain the electronic structure 
consistent with the local CEF state, we find that the 
Slater-Koster parameters for f-p hybridization should be 
limited in a certain range. A typical result is found for 
(fpn) ps and (fpcr) ~ 1 eV, which is consistent with the 
condition for the appearance of the octupole ordering. 

The organization of this paper is as follows. In Sec. II, 
in order to make this paper self-contained, we briefly re- 
view the relativistic band-structure calculations for ac- 
tinide dioxides. It is meaningful to define the problems 
included in the band-structure calculations. In Sec. Ill, 
we explain a way to construct the f-p model in the tight- 
binding approximation. Then, we determine the param- 
eters of the model, except for (fpcr) and (fpn) , from the 
comparison with the experimental and band-structure 
calculation results. In Sec. IV, we depict the energy 
band structure of the f-p model by changing the values 
of f-p hybridization. We deduce the reasonable regions 
for (fpcr) and (/p7r). In Sec. V, we discuss some future 
problems concerning the electronic structure of actinide 
dioxides. Finally, we summarize this paper. Throughout 
this paper, we use such units as h=kB=^- 



II. BRIEF REVIEW OF BAND-STRUCTURE 
CALCULATIONS FOR ACTINIDE DIOXIDES 

Let us briefly review the band-structure calculation re- 
sults in order to clarify the problem in the understanding 
of electronic structure of actinide dioxides. As for details, 
readers should consult Ref. [2ll 

We have performed the calculations by using the rela- 
tivistic linear augmented-plane-wave (RLAPW) method. 
We assume that all 5/ electrons are itinerant and perform 
the calculations in the paramagnetic phase. Note that 
we should take into account relativity even in the cal- 
culations for solid state physics because of large atomic 
numbers of the constituent atoms. The spatial shape of 
the one-electron potential is determined in the muffin- 
tin approximation. We use the exchange and correlation 
potential in a local density approximation (LDA). The 
self-consistent calculation is carried out for the experi- 
mental lattice constant for actinide dioxides. 

In Fig. 1, we show a typical result for NpC>2 along the 
symmetry axes in the Brillouin zone. In the energy band 
structure in the vicinity of Ep, there always occurs a 
hybridization between actinide 5/ and oxygen 2p states 
for actinide dioxides. The lowest six bands originate from 
the oxygen 2p states and are fully occupied and the width 
of oxygen 2p band is about 4.76 eV. Narrow bands lying 
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FIG. 1: Energy band structure of Np02 obtained by the self- 
consistent RLAPW method. Note that we pick up only 5/ 
and 2p bands around Ef, which indicate the position of the 
Fermi level. 



in the region 4.5-7.5 eV are the 5/ bands which are split 
into two subbands by the spin-orbit interaction. The 
spin-orbit splitting in the 5/ states is estimated as 0.95 
eV, which is consistent with that for isolated neutral Np 
atom. Note that in the LDA calculation, we find the 
metallic state for Np02, not the insulating state. This 
point will be discussed later. 

Here we remark that Tj doublet and Tg quartet levels 
appear around Ep at the T point. It should be noted that 
the T7 level is lower than the Ts in our band-structure 
calculations. However, from the CEF analysis on the ba- 
sis of the j-j coupling scheme, becomes lower than 
in actinide dioxides. When we accommodate 5/ elec- 
trons in r 8 orbitals, we obtain T 5 triplet for n=2, T s 
quartet for n=3, and Ti singlet for n=4, as experimen- 
tally found in the CEF ground states of UO^, NpO^, 
and Pu0 2 23 ' 24 . Note here that n denotes the number of 
local 5f electrons. 

In order to resolve the problems, it is necessary to im- 
prove the method to include the effect of CEF potentials 
beyond the simple estimation of the Madelung potential 
energy. However, it is a difficult task to perform such 
improvement concerning the formulation of the band- 
structure calculation. Thus, in this paper, we choose an 
alternative method to exploit the tight-binding f-p model 
for the purpose to understand the role of f-p hybridiza- 
tion for the change of CEF states in the tight-binding 
model. By changing the parameters in the f-p model, 
we attempt to clarify the key quantities which character- 
ize the electronic structure of actinide dioxides. 
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FIG. 2: Crystal structure of AnC>2. Large solid circles denote 
actinide ions of which positions are given by i, i + ai, i + a,2, 
and i + a-j. Small open circles indicate oxygen ions. Note 
that two oxygen ions, Ol(i) and 02(i), are included in the 
unit cell containing one actinide ion specified by i. 



III. TIGHT-BINDING APPROXIMATION 

A. Crystal structure and unit cell 

Before proceeding to the construction of a tight- 
binding model for actinide dioxides with the fluorite 
structure, first let us define the unit cell including one 
actinide ion and two oxygen ions, as shown in the Fig. 2. 
The basis vectors of the fee lattice are given by ai = 
(a/2, a/2,0), a 2 = (0, a/2, a/2), and a 3 = (a/2,0, a/2), 
where a is the lattice constant. Thus, in Fig. 2, positions 
of adjacent four actinide ions are given by i, i + a\, i + a 2 , 
and i + a 3 , where i denotes the position vector for one 
actinide ion. 

The positions of eight nearest-neighbor oxygen ions are 
given by bi = (a/4, a/4, a/4), b 2 = (-a/4, a/4, a/4), 
6 3 = (a/4, -a/4, a/4), 64 = (-a/4, -a/4, a/4), b 5 = 
(-a/4, -a/4, -a/4), 6 6 = (-a/4, a/4, -a/4), b 7 = 
(a/4, -a/4, -a/4), and b s = (a/4, a/4, -a/4). Note that 
the two oxygens, 01 and 02, in the same unit cell are 
specified by 61 for 01 and 65 for 02, respectively. 



where fima- is the annihilation operator of / electron at 
site i with spin a in the orbital specified by m. Note that 
m is the z-component m of angular momentum I = 3. We 
note also that the spin-orbit coupling is not included at 
this stage. 

Since the fluorite structure belongs to O^ point group, 
B m ,m' is given by using a couple of CEF parameters B® 
and .Bg for angular momentum I = 3 a o 25 ' 26 



-63,3 — 


B-3,-3 


= 18052 + 1805°. 


-82,2 = 


B-2,-2 


= -42052 ~ 10805°. 


Bi,i = 


B-1,-1 


= 605° + 27005° 


Bo,o = 


36052 - 


- 36005g 


B3-1 - 


= -8-3,1 


= 60/15(^2 - 215°) 


B2-2 '- 


= 30052 


+ 75605°, 



Note the relation of B m rn i = B m i, m . 

After performing the diagonalization of Hqef, we ob- 
tain three kinds of CEF states: T2 singlet (xyz), T4 triplet 
(x(5x 2 - 3r 2 ), y(5y 2 - 3r 2 ), z(5z 2 - 3r 2 )), and T 5 triplet 
(x(y 2 — z 2 ), y(z 2 — x 2 ), z(x 2 — y 2 )). The corresponding 
CEF energies are given by E(T 2 ) = -720(5° + 12.Bg), 
E(T A ) = 360(52- 10-Bg), and E(T 5 ) = -120(52-545g). 
Note that these seven states are elements of cubic har- 
monics for I = 3. In the traditional notation, we ex- 
press CEF parameters B\ and B% as B% = Wx/F(A) 
and Bl = W(l - |z|)/F(6)with F(4)=15 and F(6)=180 
for angular momentum t = 3iSS Note that x specifies the 
CEF scheme for Oh point group, while W determines an 
energy scale for the CEF potential. 

E(T 2 ) = -4SW[x + (1 - \x\)], 
£(r 4 )=4W[6or-5(l-|ar|)], (3) 
E(r s ) = -AW[2x - 9(1 - 

The value of W and x will be discussed later. 

Since we will construct the model in the cubic system, 
it seems to be natural to use these cubic harmonics as 
/-electron basis function. Thus, in the following, we de- 
fine fi as the index to distinguish the orbitals of cubic 
harmonics. Note that \i takes 1 ~ 7 and the definitions 
are as follows: 1: xyz, 2:x(5x 2 — 3r 2 ), 3:y(5y 2 — 3r 2 ), 4: 
z(5z 2 - 3r 2 ), 5:x(y 2 - z 2 ), 6:y(z 2 - x 2 ), and 7:z(x 2 - y 2 ). 
The corresponding energy E^ is given by the above equa- 
tions. 



B. CEF state 

Now we define the basis of / electrons when we consider 
the electronic model for actinide dioxides with the fluorite 
structure. For the purpose, we solve the problem of one 
/ electron in the CEF potential. The CEF Hamiltonian 
is written as 

5cef = B m>m 'f im(T fi ma , (1) 



C. Hamiltonian 

The Hamiltonian is given by 

B = Hf + Hf p + Hp, (4) 

where Hf and H p denote /- and p-electron part, respec- 
tively, while Hfp indicates f-p hybridization term. In the 
following, we explain the construction of each term. 
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1. f -electron term 



The /-electron part is given by 



and 



*t = E 



J7 



(Ef + E^S^f. fkn'a- 



(5) 



where /fc^o- is the annihilation operator of / electron with 
spin cr in the orbital fi, sj/^i is the /-electron dispersion 
due to the hopping between nearest neighbor actinide 
ions, Ef is the /-electron level, E^ denotes the CEF po- 
tential energy of /i orbital, A is the spin-orbit interaction, 
and £ is the spin-orbit matrix element. 

Concerning the expression of the spin-orbit coupling, 
it is necessary to step back to the basis of the spherical 
harmonics. On the basis labelled by m, the spin-orbit 
interaction C, m ,a,m',a' is expressed as 

Cm, ±1/2, m, ±1/2 = ±m/2, 

Cm±l/2, T l/2,m,±l/2 = \/l2 TTl{jTl i l)/2, (6) 

and zero for the other cases. By transforming the basis 
from m to we obtain Cn.o-.fj,' .a 1 in eq. ©. 

The /-electron dispersion in eq. ([5]) is expressed as 



Jf = 



E« 



ik-a±ff 



(7) 



where a denotes the vectors connecting twelve nearest 
neighbor sites of the fee lattice and vj , (a) indicates the 
/-electron hopping amplitude between fi and // orbitals 
along the direction of a. Here we note that a runs among 
±ai, ±a 2 , ±a 3 , ±(a 2 - 03), ±(03 - o{), and ±(a x - 
02). The hopping integral tjf,(a) is expressed by using 
the Slater-Koster tabl o 27 i 28 Here we consider only the 
/-electron hopping through a bond (ffcr) . 



2. f-p hybridization term 



The f-p hybridization term is written as 



H fp= E E[ y £i4^-+ h - c - 

fc,/i,f,cr j — 1,2 



(8) 



where Pjk^a is the annihilation operator of p electron 
with spin a in the orbital v of j-th oxygen and j denotes 
the label of oxygen ions in the unit cell, as shown in Fig. 1. 
Note that v runs among x, y, and z which correspond to 
p x , p y , and p z orbitals, respectively. The hybridizations 
and V^ 2 ' are, respectively, written as 



(i) 



ikai 



-ife-02 



-ik-a-i 



(9) 



he 



ik-as 



(10) 



where tyf/,(b) denotes the hopping amplitude between / 
and p orbitals along b direction. Here we note that b 
runs among 61 ~ 6§- The hopping integral t^(b) is 
represented in terms of (fpa) and (fpir) by using the 
Slater-Koster table. 2 -^ 



p- electron term 



The p-electron part is expressed as 



H P = E El £ fe£' + E P 5 io 5 uu']p\ kvrT Pjhu'a, (11) 



(ii) 

where s kllv , is the p-electron dispersion, i and j denote 
the label of oxygen ions in the unit cell, as shown in 
Fig.l, and E p is the p-electron level. Note that we take 
into account nearest neighbor and next nearest neighbor 
hoppings for p electrons. We also note that the relations 
of e 



(22) _ (11) 



kvv' 



4 J' and £ 



(21) _ J12)* 



kvv' 



"kvv' 



The diagonal part is given by 



r (n) 



2t W '(ii) cos(fc x /2 + ky/2) 
2t vv >{a,2 - a 3 ) cos(fc x /2 - ky/2) 
2t vu >(a 2 )cos(k y /2 + k z /2) 
2W'( a 3 - a i) cos(k y /2 - k z /2 
21 ' 



,/(a 3 ) cos(fc 2 /2 + k x /2) 



2t vv i{a,\ — 



a 2 )cos(fc z /2-fc x /2), (12) 



where the hopping amplitudes are given by 



tuv<{ai) — 

Uu' (02-03) = 

tuv'fa) = 

U„>(a 3 - ai) = 

t V v>(a 3 ) = 

t vv i(a\ - a 2 ) = 



P+ P- 
P- P+ 
(ppn)' 

P+ -p- 
-p- P+ 

(ppw)' 

(ppn)' 
p+ p- 
p- p+ 

(ppK)' 
P+ -p- 
-p- P+ 

P+ p_ 
(ppir)' 
P- p+ 

P+ -p- 
(ppw)' 
-p- p+ 



(13) 
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Here p± — [(ppcr)' ± (pp7r)']/2, where (jrpcr)' and (ppTt)' 
denote the Slater-Koster integral of p electron among 
next-nearest neighbor oxygen sites. 

As for the off-diagonal parts, we obtain 

e W = (ppa)[e ik < ai+a ^ + e ife ' a2 ] 
+ {ppTr)[e lk - {ai+a ^ + e ife '° 3 ] 
+ (p P 7r)[e 4fc - (a2+a3) +e u, -° I ], (14) 

egj = (^)[e 4fe -( ai+a3 )+e jfe - a2 ] 

+ (ppa)[e jfc - (ai+a2) +e jfe - a3 ] 

+ (pp7r)[e ifc - (aa+a3) +e*-° l ], (15) 

and 

4'zz = (pp7r)[e ifc -( ai+as ) +e 4fe - a2 ] 

+ (ppTT)[e tk < ai+a2) +e ife '° 3 ] 

+ (pp(i)[e 4fe -( a2+a3 ) +e a "° 1 ]. (16) 

Other off-diagonal components are all zeros. 



the ratio of the distances of nearest and next nearest 
neighbor sites, we set r\ — (I/a/2) 7 ^ 2 ~ 0.3j^ As for 
(ppcr) and (ppTr), we determine them as (ppa)-OA eV 
and (ppTr)=— 0.4 eV, after several trials to reproduce the 
structure of the wide p bands in the relativistic band 
structure calculations. 

Concerning (ffcr), it is related with the bandwidth 
W of f electrons in the j=5/2 states on the fee lattice. 
In the limit of infinite A, we have obtained W as W = 
(3/56)(50 + 2VT45)(//cr)w4.0(//cr).^Note that for the 
case of finite A, the width of J =5/2 bands is deviated 
from W, but when A is large enough as in actual actinidc 
compounds, the bandwidth is found to be almost equal 
to W. From the comparison with the relativistic band- 
structure calculation results, the width of j =5/2 bands 
is 0.5~0.7 eV, suggesting that (ffcr) is in the order of 0.1 
eV. Then, we set (ffa)=0.1 eV in the present model. 

In the following calculations, due to the diagonaliza- 
tion of the Hamiltonian, we depict the tight-binding 
bands by changing (/per) and (fpn), which are believed 
to be key parameters to understand the electronic struc- 
ture of actinide dioxides. 



D. Parameters of the model 

The tight-binding Hamiltonian includes many param- 
eters. Here we try to fix some of them from the experi- 
mental and band-structure calculations results. 

(i) CEF parameters. It should be noted that it is pos- 
sible to reproduce the CEF states of actinide dioxides, 
when we accommodate plural numbers of / electrons in 
the level scheme in which Tg is lower than Ft. As al- 
ready mentioned in Sec. II, we obtain triplet for n=2, 
Tg quartet for n=3, and Ti singlet for n=4, as experimen- 
tally found in the CEF ground states of UO^, NpO^, 
and PuG ^ 23 ' 24 . Thus, in the present tight-binding model, 
we set W — —0.01 eV and x = 0.7 in order to reproduce 
that T$ quartet is the ground state and T7 is the excited 
state with the excitation energy of about 0.2 eV. 

(ii) Spin-orbit coupling. From the relativistic band- 
structure calculation for actinide atom, the splitting en- 
ergy between j=5/2 and j—7/2 states has been found 
to be about 1 eV. Since the splitting energy is given as 
(7/2)A with the use of spin-orbit coupling A, we fix it as 
A = 0.3eV. 

(Hi) f- and p-electron levels. In this paper, the /- 
electron level Ef is set as the origin of energy, leading to 
Ef = 0. On the other hand, the p-electron level E p is 
considered to be E p — —4 eV from the comparison of the 
relativistic band-structure calculation results. 21 

(iv) Slater-Koster integrals. In the model, we use seven 
Slater-Koster integrals as {ffcr), (/per), (fpn), (ppcr), 
(ppir), (ppcr) 1 , and (ppir)'. Among them, concerning 
the p-electron hoppings, we introduce the ratio r\ be- 
tween nearest and next nearest neighbor hopping ampli- 
tudes, given by rj = (ppcr) / (ppcr)' = (ppir) / (ppir)' . From 



IV. RESULTS 

Now we show our results of the diagonalization of the 
tight-binding model. Note that in the following figures 
of the band structure, "0" in the vertical axis indicates 
the origin of the energy, not the Fermi level Ep. If it is 
necessary to draw the line of E-p , we set it from the condi- 
tion of (n)=3 for tetravalent Np ion in Np02, where (n) 
denotes the average number of / electrons per actinidc 
ion. In the present paper, we do not take care of the 
difference in actinide ions. 

First we consider the case in which the f-p hybridiza- 
tion is simply suppressed. In Fig. 3(a), we show the tight- 
binding bands for (fpa)=(fpir)=0 along the lines in the 
first Brillouin zone. We obtain the / and p bands which 
are not hybridized with each other and / bands split into 
j=5/2 and j—7/2. Note that r§ becomes lower than 
at the r point due to the effect of local CEF potentials. 
We observe some degeneracy in p bands which will be 
lifted by f-p hybridization. 

In our first impression, in spite of the simple suppres- 
sion of the f-p hybridization, the overall structure of / 
and p bands seems to be similar to that of the relativistic 
band-structure calculations in Fig. 1. However, some sig- 
nificant difference is found in the p-band structure. For 
instance, we find the level crossing in the p-band struc- 
ture of Fig. 3(a) between the L and T points, but we do 
not observe such behavior in Fig. 1. Such difference orig- 
inates from the simplification to consider only actinidc 
5/ and oxygen 2p electrons. The difference in the p-band 
structure is not further discussed in this paper. 

Next we include the f-p hybridization as (fpa)=l eV 
and (/p7r)=0.1 eV in Fig. 3(b). Due to the effect of f-p 
hybridization, we find additional dispersion in / and p 
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(a) 2.0 




X w 



X U,K 




X w 



X U,K 



FIG. 3: Energy band structure obtained by the tight-binding 
model for (a) (fpa)=(fpw)=0 and (b) (fpa)=1.0 eV and 
(/ p7r )=0.1 eV. 



FIG. 4: Energy band structure obtained by the tight-binding 
model for (a) (fpa)=1.0 eV and (/p7r)=— 0.1 eV and (b) 
(/p<r)=1.0 eV and (/prr)=-0.6 eV. 



bands. In particular, the p-band structure becomes simi- 
lar to that in the relativistic band-structure calculations. 
In this case, we still observe that Tg is lower than Ty at 
the r point. 

Let us now consider the cases of negative (fpn) by 
keeping the value of (fpa)-l eV. In Figs. 4(a) and 4(b), 
we show the results for (/p7r)=— 0.1 eV and —0.6 eV, 
respectively. For (/p7r)=— 0.1 eV, we do not find sig- 
nificant difference in the band structure from the case 
of (/p7r)=0.1 eV. However, for {fp^)—— 0.6 eV, we find 
that T*i is lower than T$ at the T point. Regarding the 
CEF states at the T point, the f-p model with (fpa)=l 
eV and (fpn)=— 0.6 eV seems to reproduce the relativis- 
tic band-structure calculation results. Note that in the 
p-band structure, we find the level crossing of two low- 



energy bands along the line between W and L points, 
which has not been observed in the band-structure calcu- 
lation. However, as mentioned above, we do not further 
pursue the difference in the p-band structure. 

Here we turn our attention to the /-electron states at 
the r point. In the relativistic band-structure calcula- 
tions for NpC>2r^ we have already pointed out that the 
level becomes lower than that of Tg, in sharp con- 
trast to the local CEF state in the j-j coupling scheme 
expected from the experimental results. This is due to 
the fact that the CEF potential is not included satisfac- 
torily in the relativistic band-structure calculation. On 
the other hand, the CEF potential is included in the 
tight-binding model within the point charge approxima- 
tion and the change of the level scheme at the T point 
can be explained by the f-p hybridization. When we do 
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FIG. 5: Energy difference A between IV and Ts states at the 
r point in the j=5/2 bands as a function of (fpn) for (fpa)=0 
eV (solid curve), (fpa)=l eV (broken curve), and (fpa) = — l 
eV (dotted curve). A positive A denotes that the energy of 
IV is larger than that of Ts . 



not consider the f-p hybridization, we find that Tg level 
becomes lower than that of T7, but with the increase of 
the effect of f-p hybridization, the order of the level at 
the r point is converted. Namely, the order of T7 and 
Tg levels is determined by the competition between the 
CEF potential and the f-p hybridization. In this sense, 
the CEF potential is not included satisfactorily in com- 
parison with the f-p hybridization in the band-structure 
calculation. 

In the fluorite crystal structure of AnC>2, actinide ion 
is surrounded by eight oxygen ions in the [111] and other 
equivalent directions. Thus, the r 7 orbital is penalized 
from the viewpoint of electrostatic energy, since its wave- 
function is elongated along the [111] directions. However, 
the wavefunctions of two r 8 orbitals are expanded in the 
directions of axes. Namely, it is qualitatively understood 
that Tg level is lower than T7 one in the actinide dioxides. 

From the viewpoint of the overlap integral between ac- 
tinide 5/ and oxygen 2p electrons, we expect that the hy- 
bridization of r 7 orbital is larger than that of Tg. Thus, 
due to the effect of f-p hybridization, the r 7 level be- 
comes lower than r 8 , even if the local CEF ground state 
is Tg. When the effect of f-p hybridization is relatively 
larger than that of the CEF potential, it is possible to 
observe that r 7 is lower than Tg, as actually found in the 
relativistic band-structure calculation results. We em- 
phasize that it is one of the key points concerning the 
f-p hybridization to understand the electronic structure 
of actinide dioxides. 

In Fig. 5, we depict the energy difference A between 
the Tg and r 7 states at the T point as functions of (/ptt) 
for several values of (/per). Note that A is positive when 
Tg is lower than r 7 . For (fpa)=0, we find that A is 
positive in the region of |(/p7r)| < 1.2 eV. When we 
change the value of (/per), A is found to be maximum 
at (fpTr) = —(/per) due to the effective disappearance of 



the f-p hybridization between actinide r 7 and oxygen 2p 
electrons. 

Readers may consider that the absolute value of {fpir) 
should not be so small only for the purpose to keep the 
order of the local CEF states. However, if we increase 
the absolute value of (fpn) for (fpa) — 1 eV, we should 
remark that the /- and p-electron bands are significantly 
changed from those in the relativistic band-structure cal- 
culation results. Thus, from the viewpoints of the local 
CEF states and the comparison with the band-structure 
calculations, the reasonable parameters are found in the 
case of small |(/j>7r)| for (fpa) = 1 eV. 



V. DISCUSSION AND SUMMARY 

In this paper, we have analyzed the tight-binding 
model for An02 in comparison with the local CEF states 
and the result of the relativistic band-structure calcula- 
tions. We have concluded that \(fpir)\ should be small 
for the case of {fpa) = 1 eV in our tight-binding model 
in order to keep the CEF levels at the T point. We have 
also emphasized that such a condition coincides with that 
for the octupole ordering on the basis of the f-p models 
Namely, the condition to keep the local Tg ground state 
is consistent with the appearance of the ordering of mag- 
netic octupole which is composed of complex spin and 
orbital degrees of freedom. 

Here we provide a comment on the local CEF state in 
the band-structure calculations. As long as we perform 
the band-structure calculations with in the LDA, it is 
found that the r 7 state becomes lower than the Tg at 
the r point, in contrast to the local CEF state expected 
from the experiment. In this paper, we have proposed the 
scenario to control the effect of f-p hybridization on the 
CEF state, but it should be remarked that in the LDA 
calculation, we could not obtain insulating state corre- 
sponding to the multipole ordering for Np02 In order 
to improve this point, we need to consider the effect of the 
Coulomb interactions, but it is a serious problem. One 
way for this problem is to employ the LDA+J7 method. 
In fact, it has been reported that we the ordered state 
with octupole and higher multipoles can be reproduced, 19 
suggesting that the r 8 state is lower than T7 in the elec- 
tronic structure. The effective inclusion of the Coulomb 
interaction is an alternative scenario to understand the 
CEF state consistent with the experiments. 

Although we have not discussed the difference in elec- 
tronic structure due to the change of actinide ions in this 
paper, it is naively expected that the difference between 
Ef and E p becomes small in the order of Th, U, Np, 
Pu, Am, and Cm from the chemical trends in actinide 
ions and the previous band-structure calculations. On 
the other hand, the change of f-p hybridization among 
actinide dioxides may play more important role to ex- 
plain the effect of the difference in actinide ions. It is an 
interesting future problem to clarify the key issue to un- 
derstand the difference in electronic structure of actinide 
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dioxides. 

In summary, we have constructed the f-p model in 
the tight-binding approximation. We have determined 
the parameters by the experimental results and the rel- 
ativistic band-structure calculations. It has been con- 
cluded that the absolute value of {fpif) should be small 
for (/pcr)=l eV in order to reproduce simultaneously the 
local CEF states and the band-structure calculation re- 
sults. The small value of |(/f>7r)| is consistent with the 
condition to obtain the octupole ordering in the previous 
analysis of the f-p model. We believe that the present 
tight-binding model will useful to extract the essential 
point of the electronic structure of actinide dioxides from 



the complicated band-structure calculation results. 
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